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1. Motivation and objectives 

One of the most exciting recent developments in the field of large eddy simulation 
(LES) is the dynamic subgrid-scale model (Germano et al 1991). The dynamic 
model concept is a general procedure for evaluating model constants by sampling 
a band of the smallest scales actually resolved in the simulation. To date, the 
procedure has been used primarily in conjunction with the Smagorinsky (1963) 
model. The dynamic procedure has the advantage that the value of the model 
constant need not be specified a priori, but rather is calculated as a function of 
space and time as the simulation progresses. This feature makes the dynamic model 
especially attractive for flows in complex geometries where it is difficult or impossible 
to calibrate model constants. 

The dynamic model has been highly successful in benchmark tests involving ho- 
mogeneous and channel flows (c./. Germano et al. 1991, Moin et al. 1991, Cabot 
and Moin 1991). Having demonstrated the potential of the dynamic model in these 
simple flows, the overall direction of the LES effort at CTR has shifted toward an 
evaluation of the model in more complex situations. The current test cases are basic 
engineering-type flows for which Reynolds averaged approaches have been unable 
to model the turbulence to within engineering accuracy. Flows currently under in- 
vestigation include a backward- facing step, wake behind a circular cylinder, airfoil 
at high angles of attack, separated flow in a diffuser, and boundary layer over a 
concave surface. Preliminary results from the backward-facing step (Akselvoll and 
Moin 1993) and cylinder wake simulations are encouraging. Progress toward the 
airfoil simulations are discussed by Choi and by Jansen, while preliminary diffuser 
simulations are discussed by Kaltenbach (all in this volume). The present paper 
discusses progress on the LES of a boundary layer on a concave surface. 

Although the geometry of a concave wall is not very complex, the boundary layer 
that develops on its surface is difficult to model due to the presence of streamwise 
Taylor- Gortler vortices. These vortices arise as a result of a centrifugal instability 
associated with the convex curvature. The vortices are roughly 1/3 of a boundary 
layer thickness in diameter, alternate in sense of rotation, and are strong enough to 
induce significant changes in the boundary layer statistics. Owing to their stream- 
wise orientation and alternate sign, the Taylor- Gortler vortices induce alternating 
bands of flow toward and away from the wall. The induced upwash and down- 
wash motions serve as effective agents to transport streamwise momentum normal 
to the wall, thereby increasing the skin friction. As evidenced by the 1980 AFOSR- 
Stanford conference on complex turbulent flows, Reynolds averaged models perform 
poorly for concave curvature since the Taylor- Gortler vortices are not resolved in 
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these calculations. Historically the ad hoc corrections for the effects of curvature 
have been unsatisfactory. The objective of this work is to investigate the effective- 
ness of large eddy simulation and the dynamic subgrid-scale model for this flow. 

The simulation targets the experimental data of Barlow and Johnston (1988). 
This experiment is an ideal test case since a rather complete set of velocity statistics 
are available for several streamwise stations. 

2. Accomplishments 

A preliminary large eddy simulation of a boundary layer along a concave surface 
has been performed. The geometry and flow conditions were close to those studied 
in the laboratory by Barlow and Johnston (1988). A limited comparison with the 
experimental data has been undertaken. These items are discussed in more detail 
below. 

2 A Numerical procedure 

The computer code for this work is an adaptation of the code written by Choi and 
Moin (1993). Boundary conditions have been generalized and the dynamic subgrid- 
scale model has been added (see Kaltenbach, this volume). The incompressible 
Navier-Stokes equations are solved in a coordinate space where two directions are 
curvilinear and the third (spanwise) direction is Cartesian. Spatial derivatives are 
approximated with second-order finite differences on a staggered mesh. A fully- 
implicit fractional step algorithm is used for the time advancement. 



FIGURE 1. Computational domain. All dimensions are referred to the boundary 
layer thickness measured at the location where the curvature begins (6o). The radius 
of curvature is R = lS.l^o- 
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Figure 1 shows the geometry for the present simulation. This geometry is simi- 
lar to that used by Barlow and Johnston (1988) in their experimental study. The 
boundary layer is allowed to develop along a flat entry section approximately 106o 
in length, where So is the boundary layer thickness measured at the location where 
the curvature begins. At this point the boundary layer encounters a constant ra- 
dius of curvature bend that turns the flow through 75°. The ratio of the boundary 
layer thickness to the radius of curvature, Sq/R, is 0.055. The domain extends 6 0 
in the spanwise direction and 36q in the wall-normal direction. According to the 
measurements of Barlow and Johnston (1988), the computational domain is just 
wide enough to enclose one pair of streamwise Taylor- G or tier vortices. Periodic 
boundary conditions are applied in the spanwise direction while a no-stress condi- 
tion is applied at the upper boundary. No-slip conditions are applied at the wall. 
Turbulent boundary layer data from an independent simulation is supplied at the 
inflow boundary (see Section 2.2 below). A convective boundary condition is ap- 
plied at the outflow. The computational grid contains 178 x 40 x 32 points in the 
streamwise, normal, and spanwise directions respectively. The mesh is stretched in 
the wall-normal direction and uniform in the other two. The grid spacings, based 
on wall units at the location where the curve begins, are Aa:"*" = 98, Ay^ in = 1, 
and A z+ = 16. 

The flow conditions match those in the experiment. The momentum thickness 
Reynolds number at the start of curvature is 1300. The experiment was conducted 
in water and is therefore incompressible. 

2.2 Inflow boundary data 

A spatially- evolving simulation such as this one requires the specification of in- 
stantaneous turbulent data at the inflow boundary. Although some level of approx- 
imation must be made, accurate inflow data is desired to insure minimal transients 
and realistic turbulence within a short distance downstream of the inlet. Fairly real- 
istic instantaneous inflow data is generated via an auxiliary large eddy simulation of 
a parallel flow boundary layer. The grid used for the inflow simulation is a truncated 
version of that used in the main simulation. It extends only one boundary layer 
thickness in the wall-normal direction and 5<5o in the streamwise direction. Periodic 
boundary conditions are applied in the streamwise and spanwise directions, while 
a no-stress boundary condition is applied at the upper boundary. The inflow sim- 
ulation is run in parallel with the main simulation in a time-synchronous fashion. 
At each time step, the velocity field is extracted from the central y — z plane in the 
inflow simulation. This data is used as the inflow boundary conditions. In practice, 
the inflow simulation can be either run at the same time as the main simulation 
or run ahead of time and the inflow data stored on disk. The inflow simulation 
increases the overall cost of the main simulation by less than 4%. 

2.3 Preliminary results 

Although the simulation was patterned after the experiment of Barlow and John- 
ston, there is one important difference in geometry between the two. The experiment 
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of Barlow and Johnston was conducted in a duct where the boundary layer thick- 
ness, 6 0 was about 1/3 of the duct width. The duct was of variable cross-section 
with the width tailored to minimize the streamwise pressure gradient on the concave 
wall. The simulation can be viewed as taking place in a constant-width duct where 
the boundary layer on the convex wall is not present. 

Since the simulation was performed in a constant width duct, the boundary layer 
will experience a streamwise pressure gradient in the vicinity of the onset of cur- 
vature. The reason for this is that a normal pressure gradient is required in the 
curved section to balance the centrifugal force associated with the streamline curva- 
ture. The normal pressure gradient requires the pressure to be higher at the concave 
wall and lower at the convex wall. The development of the normal pressure gradi- 
ent in the region where the curvature begins induces an adverse pressure gradient 
along the concave wall. This effect is shown in Figure 2. The pressure is seen to 
increase abruptly along the concave wall in the vicinity of the onset of curvature. 
The gradual drop in pressure over the entire length is due to acceleration of the 
free-stream by the thickening of the boundary layer. Also shown in Figure 2 is the 
concave wall pressure estimated from the experimentally measured velocity profiles. 
By virtue of contouring the convex wall, Barlow and Johnston were able to achieve 
a nearly constant pressure distribution. 
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FIGURE 2. Pressure distribution along the concave wall. : simulation; ■ : 

experiment. The curve begins at x = 0. 


As discussed below, an inviscid analysis has been performed to determine the 
location of the streamline in the Barlow and Johnston experiment that lies near the 
center of the duct. In future simulations this streamline will be used as the upper 
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boundary where no normal velocity and no stress conditions will be applied. For 
the present, a rough comparison with the data of Barlow and Johnston can still be 
made provided that the region affected by the pressure gradient is excluded. 
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Spanwise distance, z (cm) 

FIGURE 3. Instantaneous velocity vectors in the cross-flow plane at the 60° 
(r = 142 cm) station. The label 6 indicates the boundary layer thickness. 

Instantaneous velocity vectors in the cross-flow plane at the 60° (r = 142 cm) 
station are shown in Figure 3. A pair of Taylor-Gortler vortices is evident in the 
lower 1/3 of the boundary layer. As in the experiment of Barlow and Johnston, 
the vortices develop a few boundary layer thicknesses downstream of the onset of 
curvature and are coherent from that point to the 60° station. The vortex diameter 
is about 1/3 of the boundary layer thickness at the 60° station. The vortices enhance 
turbulent mixing near the wall find increase the skin friction as a consequence. This 
effect is shown in Figure 4 where the skin friction (presented as friction velocity) 
is plotted as a function of distance along the wall. The curvature begins at x = 0. 
The skin friction initially drops sharply as the boundary layer enters the curve, but 
this effect is due to the adverse pressure gradient in that region. Following a quick 
recovery, the skin friction undergoes a monotonic growth with streamwise distance. 
After the flow has been turned through 30° (x = 71 cm), the pressure gradient is 
minimal and the simulation results are expected to differ from the experiment only 
through history effects. Indeed the simulation results are in reasonable agreement 
with the experimental data between 30° and 60°. It is also important to note that 
the skin friction is well behaved near the inflow boundary and does not exhibit a 
visible transient arising from the approximate turbulent inflow data. 
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Streamwise distance, x (cm) 


FIGURE 4. Computed skin friction. : simulation; ■ : experiment. The 

curvature starts at x = 0. 

2.4 Streamline for zero pressure gradient 

The streamwise pressure gradient discussed above must be removed before a 
detailed comparison can be made with the experimental data. An obvious way 
to do this would be to simulate the exact geometry used in the experiment. The 
drawback of this approach is that the boundary layer on the convex wall would have 
to be resolved and consequently nearly twice as many grid points would be required. 
A more economical approach is adopted here where the location of the streamline 
that lies approximately midway between the two walls of the experimental geometry 
is determined. This streamline forms the upper boundary in the simulation and 
results in a pressure distribution very close to that found in the experiment. 

The location of the streamline is found using a procedure analogous to that used 
in designing the experimental facility. An inviscid analysis is used to determine 
the pressure distribution along the concave wall of a constant- width duct. The 
shape of the streamline that forms the duct convex wall is then iteratively adjusted 
in an attempt to minimize the pressure gradient on the concave wall. Once the 
optimal geometry is determined, the displacement thickness of the boundary layer 
is estimated and the streamline position adjusted to allow for the thickening of the 
boundary layer. 

Figure 5 shows the streamline determined by the inviscid analysis. The cor- 
responding pressure distribution along the convex wall is shown in Figure 6. For 
reference, the pressure distribution for a constant-width duct is also shown in Figure 
6. It can be seen that the pressure gradient is greatly reduced, but not completely 
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Figure 5, Streamline location for minimization of the pressure gradient at the 
wall. : streamline; : constant- width reference. Dimensions are in cm. 

eliminated. The residual excursions could not be eliminated with a simple algo- 
rithm that adjusts the streamline locally in response to the pressure deviation on 
the opposite side of the duct. Nearly identical pressure excursions are also present 
in the inviscid analysis used to design the experimental facility. For this reason no 
attempts were made to further refine the streamline. The boundary layer displace- 
ment thickness distribution was estimated from the experimental measurements 
and the streamline was displaced away from the concave wall accordingly. The 
final streamline will be used as the upper computational boundary in all future 
simulations. 

3. Future plans 

Future work will focus on refining the simulations and in making detailed com- 
parisons with experimental data. The upper boundary location will be changed as 
described above in order to minimize the stream wise pressure gradient. The mesh 
spacings will be varied in order to determine the minimal resolution for which ac- 
ceptable results are obtained. The spanwise extent of the domain will be enlarged 
so that the spacing of the Taylor- Gortler vortices are not imposed directly. Ideally 
the spanwise extent should be large enough to support several pairs of vortices. 
The spanwise length will be made as large as practical given the computer resource 
constraints. Detailed comparisons will be made between the LES, experimental 
data, and a simulation run with no subgrid-scale model. The latter will be used 
to determine the influence of the subgrid-scale model in the overall accuracy of the 
simulation. 
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FIGURE 6. Wall pressure distribution for the streamline shown in Figure 5. 
using the streamline boundary; using a constant-width boundary. 
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